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Abstract 

We study the charmful three-body baryonic B decays B — * D^NN: the color-allowed modes 
W - £>(*)+ np and the color-suppressed ones B — ► D^°pp. While the D* + /D + production ratio 
is predicted to be of order 3, it is found that D°pp has a similar rate as D*°pp. It is pointed out that 
B° — ► D(D*)NN are dominated by the nucleon vector current or by vector meson intermediate 
states, whereas 5° — ► D°pp proceeds mainly via the exchange of the axial- vector intermediate state 
ai(1260). The study of the TViV invariant mass distribution in general indicates a threshold baryon 
pair production; that is, a recoil charmed meson accompanied by a low mass baryon pair except 
that the spectrum of D°pp has a hump at large pp invariant mass m p p ~ 3.0 GeV. 
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I. INTRODUCTION 



Previously CLEO has reported the first observation of color-allowed charmful baryonic 
B decays with sizable rates 0: 

B(B° -> D* + np) = (14.5jto ± 2-7) x 10" 4 , 
B(B° -> D* + pp7i-) = (6.5±i;l ± 1.0) x 10~ 4 . (1) 

This, when combining with the non-observation of the two-body baryonic i? decays such as 
Z? — > NN 0, implies the dominance of multi-body final states in decays of B mesons into 
baryons. Recently, Belle announced a similar measurement for color-suppressed baryonic 
decay decays at the level of 1CT 4 ||: 

B(B° -> D*°pp) = (1.201°;! ± 0-21) x 10~ 4 , 

B(B° -> = (1.18 ± 0.15 ± 0.16) x 10" 4 , (2) 

with 5.6cr and 12a statistical significance respectively. Roughly speaking, the D*°pp rate is 
smaller than that of D* + np by one order of magnitude. 

Another class of charmful baryonic B decays is B — > A c (£ c )iVX. The early CLEO 
measurement |4| and the new Belle H and CLEO || results show that the three-body 
charmful decay B~ — ► A c pTT~ has a magnitude larger than B° — > A c p. These modes have 
been theoretically studied in ]7J. The recent first observation of the penguin-dominated 
charmless baryonic decay B~ — > ppK~ by Belle || clearly indicates that it has a much 
larger rate than the two-body counterpart 5° — > pp. Theoretically, it has been explained in 
why some charmless three-body final states in which baryon-antibaryon pair production 
is accompanied by a meson have a rate larger than their two-body counterparts. 

Under the factorization assumption, the decay amplitude, dominated by the color-allowed 
external lU-emission, is proportional to ai(Oi)f act where 0\ is a charged current-charged 
current 4-quark operator, while the decay amplitude, governed by the factorizable color- 
suppressed internal PF-emission, is described by a2(02)fact w ith O2 being a neutral current- 
neutral current 4-quark operator. Since 5° — > D^ + np are color-allowed, while 5° — > 

£>(*)+ 

np are color-suppressed, it is naively expected that the measured ratio R = T(B — > 
D^ pp)/T(B° — > D^ + np) can be used to extract the parameter \a 2 /ai\, just as in the case 
of B — > D°n° and B — > D + n~ decays. However, there is one complication here: the 
factorizable decay amplitude of color-suppressed baryonic B decays B — > D^°pp involves 
a three-body hadronic matrix element which is basically unknown. Therefore, one needs 
to impose further theoretical assumptions in order to extract a 2 from the color-suppressed 
baryonic B decay modes. The color-favored decays B° — > D^ + np have been studied in 
[TT| . It turns out that D* + / D + production ratio is predicted to be of order 3. It is thus 



anticipated that the D*°/D° production ratio in color suppressed decays is also of order 3. 
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However, experimentally the latter is consistent with unitary. This motives us to investigate 
why D°pp has a similar rate as D*°pp. 

All B -> D&NN decays can be described in terms of the pole model; they receive 
contributions from various intermediate states: vector mesons such as p, u, axial- vector 
mesons ai(1260), /i(1285), /i(1420), and pseudoscalar mesons ir,r),r)'. It appears that the 
decay B — > D°pp is very special: it is dominated by the axial-vector meson states, whereas 
the other modes proceed mainly through the vector meson poles. This enables us to explain 
the similar rates for D°pp and D*°pp. 

This paper is organized as follows. We first discuss the color-favored modes B° — > D^ + np 
in Sec. 2 and then turn to the color-suppressed ones B° — ► D^°pp in Sec. 3. Discussions 
and conclusions are presented in Sec. 4. 



II. COLOR- ALLOWED if D^+np 

At the quark level, the color-allowed decays B° — > D + ^np proceed through the fac- 
torizable external H^-emission and jy-exchange diagrams, and the nonfactorizable internal 
H^-emission [see Fig. 1(a)], while W -> £>°W pp via the factorizable internal VT-emission, 
VF-exchange diagrams and the nonfactorizable internal jy-emission [see Fig. 2(a)]. More 
precisely, their factorizable amplitudes read 

A (B° _> D+^np) iact = ^Kr rf ^{a 1 (np|(^)|0)( J D + W|(c6)|5°) 

+ a 2 (D + ^np\{cu)\0)(0\{db)\B°)}, 
A(B° -> D°^pp) iact = ^V* d V cb {a 2 {np\ (db) \B°) (D ^ \ (cu)\0) 

+ a 2 (D^pp\(cu)\0)(0\(db)\B°)}, (3) 

where (qiq 2 ) = 9i 7/^(1 — 75)92, and ai, a 2 , which will be specified later, are some renormaliza- 
tion scale and scheme independent parameters . The second term in each decay amplitude of 
Eq. (|3|) corresponds to the ^-exchange amplitude, which is not only color but also helicity 
suppressed. Since the three-body matrix element (np\(db)\B ) is basically unknown, we will 
evaluate the color-suppressed amplitude based on the pole model approximation. 

Let us first focus on the decays B — > D + ^np. To evaluate their factorizable amplitude, 
we need to know various haronic matrix elements. The one-body and two-body mesonic 
matrix elements are given by [|l(| 



(P(g)K|0> = -i/ P g M , (Vip,e)\V^\0) = f v m v e; 
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FIG. 1. (a) Quark diagrams for B 
external PF-emission. 
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where q = p B - p, F x (0) = F (0), A 3 (0) = A (0), and 

2 m B + m v A fJ2s m B - m v 
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The two-body baryonic matrix element appearing in Eq. (|3]) can be parametrized as 



(n(p n )P(pp)\(V - A) M |0) = u n {p n )l f7{q 2 )l, + • 



npi 2\ i ■92 P ( ( 1 2 ) 



„ n Pl „2\ 
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2m 
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2m 
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75 \vp{pp), 



(6) 



with q = p p + pp. Among the six baryonic form factors, the vector form factor fz{q 2 ) 
vanishes because of conservation of vector current (CVC), while the absence of the second- 
class current implies (?2(<? 2 ) = 0. Using CVC, the vector form factors /i,2(<? 2 ) can be related 
to the electromagnetic form factors of the nucleon defined by 
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FIG. 2. (a) Quark diagrams for 5° — > D°pp and (b) the pole diagrams for the factorizable 
internal VF-emission. 



F2(q 2 ) v 



(7) 



(iV(p 1 )iV(^)|^ m |0} = ^(pi)[F 1 (g 2 ) 7M +i ^ 
Specifically (see e.g. |J), 

/i^) = ^ 2 (*)-^), (8) 

where t = q 2 . 

The experimental data of the nucleon's e.m. form factors are customarily described in 
terms of the electric and magnetic Sachs form factors G^(t) and G^(t) which are related to 
Ff and F? via 



G p nt) = pr(t) + j^pnt), G% n (t) = Fnt) + 



(9) 



A recent phenomenological fit to the experimental data of nucleon form factors has been 
carried out in JO] using the following parametrization: 



t 2 t 3 t 4 t 5 t 6 



111 m 



\GUt)\ = fe + n) (in 

where Q = A QCD ~ 300 MeV and 7 = 2 + ^ = 2.148 . We will follow (TJ to use the best 
fit values 



t 

Ql 



(10) 
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xx = 420.96 GeV 4 , x 2 
x A = -433916.61 GeV 10 , 



-10485.50 GeV , 
x 5 = 613780.15 GeV 



x 3 = 106390.97 GeV 8 
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(11) 



and 



yi = 236.69 GeV 



y 2 = -579.51 GeV 



(12) 

extracted from the neutron data under the assumption \G E \ = Note that the form 

factors given by Eq. fllCf ) do satisfy the constraint from perturbative QCD in the limit of 
large t [|TTJ. Also as stressed in [O], time-like magnetic form factors are expected to behave 
like space-like magnetic form factors, i.e. real and positive for the proton, but negative for 
the neutron. 

In terms of the nucleon magnetic and electric form factors, the weak form factors read 



fi P (t) 

f7(t) 



4ml 



-G p M (t) - G E {t) 



4m' 



-G n M {t) ~ G%(t) 



N) 



t/(4m 2 



N ■ 



G p M (t) - G p E (t) + G n M (t)-Gl(t) 



t/(4m 



1 



t/(4m 



1 



(13) 



According to perturbative QCD, the weak form factors in the large t limit have the asymp- 
totic expressions |1| 



f7(t)^G p M (t)-G n M (t), g7(t) 



lG p M (t)+G n M (t). 



(14) 



It is easily seen that this is consistent with the large t behavior of fi P given by Eq. (|13l) . 

For the axial form factor gi(t), we shall follow [|L4] to assume that it has a similar 
expression as G^it) 



97(t) 



In 



Qo 



(15) 



where the coefficient d\ is related to x\ and y\ by considering the asymptotic behavior of 
Sachs form factors G P M and G 7 ^ [see Eq. (|T0|)1 



-xi - yi- 



(16) 



As shown in ||, the induced pseudoscalar form factor g 3 corresponds to a pion pole contri- 
bution to the np axial matrix element and it has the form 



t - mi 



(17) 
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III. COLOR-SUPPRESSED B° -► D^°pp 



We next turn to the color-suppressed modes B — > D^°pp and assume that the main 
contributions arise from the factorizable internal W^-emission diagram [see Fig. 2(a)]. There 
are two corresponding pole diagrams as depicted in Fig. 2(b): one with bottom baryon poles 
Sft(| + ) and ), and the other with meson poles. The low-lying meson intermediate 

states are: 7r°, 77, 77', p°,u, and a? and possibly /1 (1285) and /i(1420). 



A. Meson pole contributions 



The meson pole contribution from Fig. 2 (b) is 



A(B° -> D°pp) M = ^V: d V cb a 2 {D»\(cu)\V){ - [(A(db)\B°) 
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(18) 



where g = jj b — p D = p p + Pp- For simplicity, we have concentrated on the low-lying poles 
and neglected those contributions from the higher axial vector meson states such as /i(1285) 
and /i(1420). As we shall see, the vector and tensor coupling constants g pNN and g% NN 
are related to the vector form factors fi P and f 2 p respectively, while g1 lNN and g^NN are 
connected to the axial-vector form factors g™ v and g% p respectively. 
After some manipulation we obtain 
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(19) 
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where we have applied the relations u p $v p = 0, (p p — pp) ■ q = and employed the form 
factors defined by 

(a^p ai )\(ub) v _ A \B°(pB)) = ^W^^l? 2 ) 



m B + m ai 



Bau„2\ £ -PB (m i „ \ T^Sp/ 2^ 



+ m ai 
V 3 Bai (q 2 )-V Ba ^q 2 ) 



(20) 



with 



2 mg+% 2 m B - my 2 n 



(21) 



and V3(0) = Vo(0). Note that the pion in the form factor F Bn is referred to the charged one, 
so are the form factors Aq P and V Bai . 

Likewise, the meson pole contribution to B° — > D*°pp reads 

A(B° -> D*°pp) M = %V: d V cb a 2 f D ,m D J\V2F^(ml f )^\ 

V 2 [ L q —m* 



9r)'NN 



q 2 — m 2 



q 2 — m 2 , 
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(22) 



B. Baryon pole contributions 

In addition to the aforementioned meson pole contributions, there also exist baryon pole 
diagrams, namely, the strong process B° — > S^^V followed by the weak decay S^*- 1 — > _D°p. 
Due to the large theoretical uncertainties with the \ state S^*, we will focus on the | + 
intermediate state and its amplitude is given by 

G F . 1 



X M p 



{fi bP ijn 2 D )[2p D ■p p + p , D ( y m J]b -m p )]7 5 
+ 9i bP {m 2 D )[2p D ■ p p - ^> D (m^ b + m p )]}v p , (23) 
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where we have applied the factorization approximation to the weak decay — > D°p. Sim- 
ilarly, for B° — > D*°pp we have 

G 1 
A(B° - D«pp) B = ^V cb V: d f D ,m D * a 2 g^ Bp ^ +pfl . )2 _ m 2, 

+ ^f bP (^i)*)[ 2 PpM - ( m s b + m p )7 M + -f^ D }}v p , (24) 

where is the polarization vector of the D*. 

The baryon pole contribution is expected to be suppressed relative to the meson pole due 
to the smallness of the strong coupling of — ► B°p |§ . 

IV. CALCULATIONS AND RESULTS 

To proceed numerical calculations we first need to know the relevant form factors, decay 
constants, strong couplings, and the parameters ai, a 2 , which will be discussed in more detail 
below. 



A. form factors and decay constants 



For the mesonic form factors of B — ► P and B —>■ V transitions we use the Melikhov-Stech 
(MS) model based on the constituent quark picture and the Neubert-Rieckert-Stech-Xu 
(NRSX) model [16] which takes the Bauer-Stech-Wirbel (BSW) model [[HJ results for the 
form factors at zero momentum transfer but makes a different ansatz for their q 2 dependence, 
namely, a dipole behavior is assumed for the form factors Ft, A , A 2 , V, motivated by heavy 
quark symmetry, and a monopole dependence for F ,Ai. 

For B — > di form factors, there are two existing calculations: one in a quark-meson 
model [17[ and the other based on the QCD sum rule []18[ . The results are quite different, 
for example, V Bai (0) computed in the quark-meson model, 1.20 , is larger than the sum-rule 
prediction, —0.23 ± 0.05 , by a factor of five. We shall see later that B° — ► D°pp is rather 
sensitive to the form factor V Bai . It turns that in order to accommodate the measurement 
of this decay, this form factor should be around 0.86 which is between the above-mentioned 
model calculations. In the present paper we shall use the quark-meson model results for the 
B — > ai form factors except that the value of V Ol (0) is replaced by 0.85 rather than 1.20. 

To compute the form factors for Fq V and Fq V , it is more natural to consider the flavor 
basis of T] q and rj s defined by 



1 

V2 



'uu + dd), 



Vs 



ss. 



(25) 
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The wave functions of the 77 and rj' are given by 

, <h\ in. \ 

(26) 




cos (p — sm 
sin 6 cos 




where <p = 9 + arctanA/2, and 6 is the 77 — 77' mixing angle in the octet-singlet basis. The 
physical form factors then have the simple expressions: 



F$ = ^= cos F*r , F B f = -±= sin . (27) 



Using i^f^O) = 0.307 and Fq V ' uU (0) = 0.254 obtained from and the mixing angle 
= 39.3° (or 9 = -15.4°) || we find Fq V (0) = 0.168 and F W (0) = 0.114 in the BSW 
model and hence the NRSX model. For other form-factor models, we shall apply the relation 
based on the isospin-quartet symmetry 

^0,1 - -^0,1 - ^0,1 { Z °J 

and Eq. (p7| ) to obtain the physical B — > 77 and B — > 77' transition form factors. For the MS 
model we obtain Fq V (0) = 0.141 and F S,? '(0) = 0.115. 

For the heavy-light baryonic form factors f 1 b and g i b , we will follow |20[ to apply 



the nonrelativistic quark model to evaluate the weak current-induced baryon form factors 
at zero recoil in the rest frame of the heavy parent baryon, where the quark model is most 



trustworthy. Following [21] we have 



fx tP (q 2 J = 1-703, gf tp (q 2 J = -0.166 (29) 

at zero recoil q 2 ^ = (m Sb —m p ) 2 . Since the calculation for the q 2 dependence of form factors is 
beyond the scope of the non-relativistic quark model, we will follow the conventional practice 
to assume a pole dominance for the form-factor q 2 behavior: 

m=m ^Z^J, s( rt = , (A) (lz£Mi)\ ,30) 

where my {ttia) is the pole mass of the vector (axial-vector) meson with the same quantum 
number as the current under consideration. 

For the decay constants we use f D = 200 MeV, f D . = 230 MeV and /„, = 205 MeV. 



B. strong couplings 

In order to compute the decay rate for B° — > D^°pp we also need to know the strong 
couplings g-xNN, g\ N , #2 N an d gi lNN and their q 2 dependence. To do this, let us con- 
sider the pole contributions to B — > D^ + np. In the pole model description, the relevant 
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intermediate states are tt , p and a 1 (1260) as shown in Fig. 1(b). The matrix element 
(np\(V — A) M |0) then reads 



(np\(V - A) M |0)p O i e = u n 



- (V2f P m p [ NN .gf N 



1 



i p m p 
q 2 — m 2 



2 — 9i 7m7s ^ 0—^75 K- 31 

g z — m~ q z — J 

Comparing this with Eq. (^) we see that the p~ meson is responsible for the vector form 
factors fi and f2, % (1260) for gi and g^ and 7r~ for the induced pseudoscalar form factor 
g 3 . More precisely, 



V2f p m p Ji v * " a/2/. 



/ p"<>p V pi/lp 

tf^V) = ^T^H? 2 ), ftrW) = f " m ' } ggV)- (32) 

V2f ai m ai 2\/2/ vr m A r 

As for the vector and tensor couplings of the u meson, a priori they are not necessarily related 
to those of the p meson. For simplicity we shall assume that gi2 N (l 2 ) = 9\% N {q 2 ) I noting 
that the p meson here is referred to the charged one. 

As for the strong rj and v[ couplings with nucleons, we shall apply the 3 Pq quark-pair 
creation model |]22|j23|] to estimate its strength relative to the pion. This model in which the 
qq pair is created from the vacuum with vacuum quantum numbers 3 Po implies 

9vpp = (Ml24)$,(35)|Ml23)$ vac (45)) (33) 
g^ pP ($ p T(124)$ 7r o(35)|$ pT (123)$ vac (45))' 1 ' 

where the $'s are the spin-flavor wave functions and the vacuum wave function has the 



expression 



$ vac = -j={uu + dd + ss)®-j=(U + IT). (34) 
Using the proton wave function 

p T = -J=[^ Xs + (12) + (13)], (35) 

with abcXs = {2a}Wc^ — aW<r — a) ^c T )/\/6, the 7T° meson wave function 

= ±={uu - dd) ® -^(U - IT), (36) 

and the 77 and rf flavor wave functions given by Eq. (p6[) , we obtain 

9rjNN 3 9n'NN 3 . 

— ■ = -cos0, — ■ = -sin0. (37) 

g-nNN 5 g-KNN 5 
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Strictly speaking, the above relations hold only at low energies. But we shall assume their 

validity at arbitrary q 2 . 

As for the strong coupling g + -^o , we use the experimental result for B~ — > A c pn~ to 

H ^ B p 

fix the absolute coupling strength of 9 A +_+-]f which is in turn related to 5 , E +_ > ^° p via the 3 Pq 
quark-pair-creation model 0. It is found that \g^ + ^o \ ~ 5. 



C. a± and a 2 



In the naive factorization approach, the parameters a\ and are given by a,\ .2 = c 2j i + 
01,2/ N c , but this does not include nonfactorizable effects which are especially important for 02- 
Phenomenologically, one can treat as free parameters and extract them from experiment. 
The experimental measurement of B — > J/ipK leads to \a,2(J / ipK)\ = 0.26 ± 0.02 This 



seems to be also supported by the study of B — ► Dn decays: Assuming no relative phase 
between a\ and a 2 , the result a 2 ~ (9(0.20 — 0.30) |24],|25j is inferred from the data of 
B° — > D^ + -k~ and B~ — > D^°n^ . However, the above value of 02 leads to too small decay 
rates for 5° — > D^°tt° when compared to recent measurements by Belle and CLEO ||26|| . 
In order to account for the observation, one needs a larger 02 (-Dtp) with a non-trivial phase 



relative to a\ |27H29|j 



Using the measurements of CLEO and Belle for B° — ■> _DM°7r [26[j, the magnitudes of 



a\ and 02 and their relative phase are extracted in [28|, as exhibited in Table I. We will use 
the values of ai^(Dir) to compute the decay rate for B — > D + np, D°pp and ai^{D*%) for 



B° -> D*+np, L>*>p. 



TABLE I. Extraction of the parameters a\ and 02 from the measured B —* D^ir rates by as- 
suming a negligible PL-exchange contribution. Note that a,2(D7r) and a 2 (-D*'7r) should be multiplied 
by a factor of (200 MeV/f D ) and (230 MeV/f D *), respectively. This table is taken from [28]. 



Model 




|a 2 (D7r)| 


a 2 (-D7r)/ai(D7r) 


|ai(_D*7r) 


|a 2 (D*7T)| 


a 2 (-D*7r)/ai(D*7r) 


NRSX 
MS 


0.85 ± 0.06 
0.88 ±0.06 


0.40 ± 0.05 
0.47 ±0.06 


(0.47 ±0.05) exp(i59°) 
(0.53 ±0.06) exp(i59°) 


0.94 ± 0.04 
0.85 ± 0.03 


0.31 ± 0.04 
0.39 ±0.05 


(0.33 ± 0.04) exp(i63°) 
(0.46 ± 0.06) cxp(i63°) 



D. Results 



In principle, the unknown parameter d,2 appearing in the form factor g\ P {q 2 ) [cf. Eq. 
(|15|)1 can be fitted to the measured central value of the branching ratio for B — > D* + np as 
it is theoretically much more clean. However, we find that the decay rate of B° — > D°pp 
is dominated by the axial-vector meson poles and hence it is rather sensitive to gi P {t) and 
hence d.2- Therefore, we instead fix it by fitting to the measured central value of E{B — > 
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D°pp) = 1.18 x 1(T 4 . We obtain d 2 = -2070 GeV 6 and -2370 GeV 6 , respectively, in NRSX 
and MS form-factor models. 

The total decay rate for the process B(pb) — > N(px) + N(p 2 ) + D(p 3 ) is computed by the 
formula 

T = - — — 5- f \A\ 2 dm^drr3,n, (38) 

(2tt) 3 32m| J K J 

where mfj = (pi +Pj) 2 with p 3 = po- To compute the branching ratios, we use the B meson 
lifetimes quoted in f3(|. 

TABLE II. Branching ratios (in units of 10 -4 ) for charmful decays 5° — ► D^ + np and 
B° — > D(*) pp calculated in the MS and NRSX form-factor models. The first (second) number 
in parentheses is the branching ratio due to the vector (axial-vector and pseudoscalar) current or 
intermediate vector (axial-vector) meson contributions. 



Decay MS NRSX Expt. g],§ 

£ U D+np 3.17 (3.04,0.12) 3.64 (3.47,0.16) 

B° -> D*+np 10.0 (9.54,0.49) 11.0 (10.2,0.76) 14.5i| 4 ± 2.7 

B° -> D°pp 1.18 (0.15,1.03) 1.17 (0.11,1.06) 1.18 ± 0.15 ± 0.16 

B° -> D*°pp 1.58 (1.42,0.17) 1.23 (1.12,0.11) 1.20±g^ ± 0.21 



The results are shown in Table II. As stated before, we fit the unknown parameter d 2 to 
the measured branching ratio of B° — > D°pp and then in turn predict other neutral baryonic 
B modes. The baryon pole contributions to D^°pp are found to be at most of order of 10~ 6 
and hence they are negligible. It is clear from Table II that the predicted rates are consistent 
with experiment. We see that S° — > D(D*)NN are dominated by the vector current or by 
vector meson intermediate states,[] whereas B° — > D°pp is dominated by the axial-vector 
intermediate state ai(1260). Note that the ratio D* + / D + is of order 3, while D*°/D° is close 
to unity. 

In Fig. 3 we show the np invariant mass distributions dB/dm n p of B — > D* + np and 
B — > D + np, where m n p is the invariant mass of the nucleon pair. Evidently, the spectrum 
peaks at m p p ~ 2.1 GeV, indicating a threshold enhancement for baryon production, that 



N As far as the axial vector contribution to the branching ratio is concerned, our result for B 



D* + np is quite different from that given in [14|, though the value of d 2 is similar. For example, 
Ba ~ 12.7 x 10 -4 is obtained in [14|, whereas it is only 0.6 x 10 -4 in our case. If g™ p is identified 
with the asymptotic form \G P M + G 1 ^ in the whole time-like region [cf. Eq.(|l4|)1, Ba in our case 
will be of order only 8 x 10 -6 , while it can be as large as 1 x 10~ 4 in [14]. 
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FIG. 3. The np invariant mass distribution dB/dm np of (a) B° — > D* + np and (b) B° — > D + np. 

is, a recoil charmed meson is accompanied by a nucleon pair with low invariant mass. This 
effect is due to the suppression of the baryonic form factors at large t. Physically, this can 
be visualized that the quark and anti-quark forming a nucleon pair are moving collinearly 
and energetically, so that the invariant mass m pp tends to be small and near threshold. 

Also shown in Fig. 4 are the pp invariant mass distributions of 5° — > D*°pp and 5° — ► 
D°pp calculated in the MS model. It is clear that the spectrum of D*°pp is similar to that of 
]j{*)+ n p. As for the differential rate of D°pp, it has a hump around m pp ~ 3.0 GeV, which 
is caused by the mass term q fl q u /m 2 ai in the propagator of the a\ meson]] We see that the 
predicted spectrum is consistent with experiment. 

V. CONCLUSIONS 

We have studied the charmful three-body baryonic B decays B — > D^NN: the color- 
allowed modes B — > D^ + np and the color-suppressed ones B u _> D^°pp. While the 
D* + /D + production ratio is predicted to be of order 3, it is found that D°pp has a similar 
rate as D*°pp. It is pointed out that B — > D(D*)NN are dominated by the nucleon vector 
current or by vector meson intermediate states, whereas B — > D°pp proceeds predominately 
via the axial- vector intermediate state ai(1260). The study of the iViV invariant mass dis- 
tribution in general indicates a threshold baryon pair production, that is, a recoil charmed 
meson accompanied by a low mass baryon pair except that the spectrum of D°pp has a hump 
at large pp invariant mass m p p ~ 3.0 GeV. The presence of a hump in the D°pp spectrum 



tThe B -> D° pp spectrum is somewhat sensitive to the model for form factors. In the NRSX 
model, the peak appearing at low pp invariant mass ~ 2 GeV is lower than the hump at m p p ~ 3 
GeV. This is inconsistent with experiment [see the data shown Fig. 4(b)]. 
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mpp (GeV) mpp (GeV) 



FIG. 4. The pp invariant mass distribution dB/dm p p of (a) B — ► D*°pp and (b) B — > D°pp. 
The experimental data for the spectrum of £ — ► D°pp are taken from [3]. 

can be tested by the improved experiment in the future. 
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